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ABSTRACT 

We investigate the anisotropy of Alfvenic turbulence in the inertial range of slow solar wind 
and in both driven and decaying reduced magnetohydrodynamic simulations. A direct com- 
parison is made by measuring the anisotropic second-order structure functions in both data 
sets. In the solar wind, the perpendicular spectral index of the magnetic field is close to —5 /3. 
In the forced simulation, it is close to —5 /3 for the velocity and —3/2 for the magnetic field. 
In the decaying simulation, it is —5/3 for both fields. The spectral index becomes steeper at 
small angles to the local magnetic field direction in all cases. We also show that when us- 
ing the global rather than local mean field, the anisotropic scaling of the simulations cannot 
always be properly measured. 
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1 INTRODUCTION 

The solar wind 



the solar wind is a turbulent plasma (see reviews by 
iGoldstein et al.ll995l ; lBruno & Carbonel2005l ; lHorburv et alj2005T) 
with a power spectru m extending over many orders of magnitude 
(e.g. IColemanl 1 19680 . Scales larger than the ion gyroradius are 
known as the inertial range, and the spectral indices at 1 AU are ob- 
served to be close to —5 /3 for the magnetic and electri c fields and 
-3/2 for velocity (e.g . Matthae us~& Goldstein! 1 19821: iBale et al. 
20051; IPodesta et al.l20071;lTessein et al.l2009l:|Podesta & Borovskv 



20100 . There is also evidence that the fluctuations are pre domi- 
nantly Alfvenic (e.g. iBelcher & Davis|[l97ll ; iHorburv et al.ll 19951 ; 
IBale et al.ll2005h . 

Solar wind turbulence is anisotropic with respect to the di- 
rection of the magnetic field. For example, the magnetic field cor- 
relation length has been shown to vary depending on the angle 
of observation with respect to the field direc tion (Crooker et al. 
19821; iMatthaeus et~ai]|l990l: basso et al.|[2005l ; lOsman & Horbury 



2007; Weygand et alj |2009). The magnetic field power and spec- 
tral indices are also observed to be anisotropic: pow er at a fixed 
scale increases with ang l e to the magnetic field dBieber et al.l 
1 19961 : iHorburvetal .11 19981 ; lOsman & Horbury||2009l) and the spec- 
tral index varies from —2 at small angles t o between —3/2 and 
—5/3 in the field perpendicular direction dHorbury et alj 1200 8t 
|Podestall2009l ; iLuo & WullioTol ; Iwicks et al]|2010l . l201ll) . These 
observations are consistent with theories of critically balanced 
magnetohydrodynamic (MHD) turbulence, for example, that of 
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iGoldreich & SridhaJ { 1995), which predicts anisotropic fluctua- 
tions (kx > fell) and the —2 and —5/3 spectral indices. 

Many simulations of plasma turbulence have been performed 
most of which have used the equations of incompressible MHD 
When a strong mean magnetic field is present, the spectral index of 
the to t al energy is closer t o —3/2 than —5/3 dMaron & Goldreicl 
| 200l |; iMuller et alj |2003| : iMuller & Grappinl 120051: iMason et 
I2OO8I ; IPerez & Boldvrevl I2OO8I : iGrappin & Mulled l201Ch . al- 
though a limit ed inertial range and the bottleneck effect 
Falkovich1ll994l) make this number hard to determine precisely 
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Beresnyak & Lazarian 2009; Beresnyak 2011). 



Anisotropy has also been measured in MHD simulations. 
Early 2D simulations showed that the turbulence develops 
wavevector an isotrop y so that the fluctuations have > fell 
dShebarinetai]|l983l) and this was later confirmed in 3D simu- 



lation s dOughton et all 1 19941 ; IMatthaeus et alj|l996l ; iMilano et all 
l200lh . The anisotropy was found to be scale-dependent, 
such that jfc|| ~ k 1 ^ dCho & VishnTacl |2000|: iMaron & Goldreichl 
l200lh. in agreement wi th the critical balance predictions 
(Goldreich & Sri dhail 1 19951) . An important point noted in these 



point not ed in thes e 
; dHorbury et aTT2 008), 



studies, and also in solar wind measurements ( 
was that the anisotropic scaling is with respect to the scale- 
dependent local mea n field and not the global mea n field. 

The the o ry of IGoldreich & SridhaJ (l995) was modified 
by iBoldvrevl d2006l) by including a phenomenon called scale- 
dependent dynamic alignment. In this theory, the velocity and mag- 
netic field fluctuations align to within a smaller angle at smaller 
scales and the perpendicular spectral index becomes —3/2. There 
is eviden ce for this scale-dep endent dynamic alignment in the so- 
lar wind dPodesta et al.ll2009l) and some driven MHD simulations 



2 C. H. K. Chen et al. 



Table 1. Mean parameter values for the 65 solar wind intervals 



Solar wind speed (v sw ) 


360 ± 10 km s 1 


Ion number density («,■) 


8.6 ±0.4 cm~ 3 


Alfven speed (v^) 


40±2kms-' 


Perpendicular ion temperature (TjjJ 


7.5 ± 0.4 eV 


Ion beta (ft) 


1.1 ±0.1 


Ion temperature anisotropy (T i± /Tj^) 


0.5 ± 0.2 


Ion gyroradius (p,) 


74 ±3 km 


Alfven ratio (r&) 


0.72 ± 0.04 



(Mason et al. 2006, 2008), although higher resolution simulations 
sugge st that the alignment saturates at small scales (Beresnyak 
l20ll . 

To date, there has not been a measurement of the spectral in- 
dex parallel to the local magnetic field in simulations. Measure- 
ments of the perpendicular spectral index in the solar wind and in 
simulations are also not always in agreement. It is important to be 
sure that the same quantities are being measured in both the solar 
wind and simulations and the subject of this paper is such a com- 
parative study. We apply a similar analysis technique to both solar 
wind data and reduced MHD (RMHD) simulations, to make a di- 
rect comparison of the anisotropic scaling. In Section[2]we present 
the solar wind analysis, in Section[3]we present the simulation anal- 
ysis, in Section|4]we compare the local and global mean field meth- 
ods and in Section[5]we present our conclusions. 



2 INERTIAL RANGE SOLAR WIND MEASUREMENTS 
2.1 Data Intervals 



In this section, we apply the multispacecraft method of IChen et al.l 
l l2010al) to obtain the power and spectral index anisotropy of iner- 
tial range turbulence in the slow solar wind at 1 AU. The technique 
is applied to 65 1-hou r intervals of data from the Cluster spacecraft 
jEscoubet etaT]|200ll) from December 2005 to April 2006, when 
the typical separation between the four spacecraft was ~ 10,000 
km. The selected intervals are from the parts of the Cluster orbit 
where the spacecraft were in the free solar wind upstream of the 
bow shock at geocentric distances of between 15 Re and 20 Re- 
They contain no evidence of ion foreshock activity: signatures typ- 
ical of the ion foreshock, such as enhanced magnetic field fluctu- 
ations and high-energy ions, are not present. The time series were 
also inspected visually to ensure that they are approximately sta- 
tionary and do not contain shocks or magnetic clouds. 

In the analysis, we use 4 s measureme nts of the magnetic field 
from the fluxgate magnetometer (FGM) tealogh et alj|200lh and 
velocit y and density m oments from the Cluster ion spectrometer 
(CIS) jReme et alj[200lh . The mean values of various parameters 
for the 65 intervals are given in TableQ] The geometric mean is used 
for the ion beta, temperature anisotropy, gyroradius and Alfven ra- 
tio. The intervals are in slow solar wind with a speed < 550 km 

s-K 

The Alfven ratio is the ratio of energy in the velocity u to 
the magnetic field in Alfven units b, and can be calculated spec- 
trally, r& = E u /E b , where E" and E b are the power spectra of 
u and b. We calculate the average in the spacecraft frequency 
range 2 x 10~ 3 Hz to 1 x 10 -2 Hz, which roughly corresponds 
to scales 36, 000 km to 180,000 km under Taylor's hypothesis 
(Tavlo3 [l938l) . While this is at larger scales than the following 
anisotropy measurements, it is in the range where noise does not 
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Figure 1. Histogram of normalised cross helicity, o~ c , for 53 of the the solar 
wind intervals in the spacecraft frequency range 2 x 10~ 3 Hz to 1 X 10 -2 
Hz. 



appear to dominate the velocity spectra. The value slightly less 
than unity that we obtain (~ 0.7) is consistent with previous mea- 
surements (e g.lMatthaeus & Goldsteirj|l982l; iMarsch & Tulll99fi 
IPodesta et al]r2007l : lBruno et alfeOOllSalem et al.ll2009h . 
We also calculate the normalised cross helicity, 



E+-E' 
E++E- 



(1) 



where E + and E~ are the power spectra of the Elsasser variables 
Tr^ = u ± b. The average value for each interval is calculated over 
the same range as the Alfven ratio. The usual convention is used: 
the Elsasser variables are defined such that positive values of <7 C 
correspond to Alfvenic propagation away from the Sun. A his- 
togram of a c (Fig.[T]( shows a range of values with a non-Gaussian 
distribution: there is a large outward population (o~ c > 0.5), a bal- 
anced population (ct c ~ 0), and a few inward intervals (cy < —0.5). 



2.2 Analysis Technique 

For each interval, pairs of points from the time series of the four 
spacecraft are used to calculate second-order structure functions 
at different angles to the local magnetic field, as described by 
IChenetalJ d2010a). The second-order structure function is defined 



8Bf(l) = (\B i (r + l)-B i (T) 



(2) 



where 6, is the fth component of the magnetic field, 1 is the sepa- 
ration vector, and the angular brackets denote an ensemble average 
over positions r. The local mean magnetic field at scale 1 is defined 

as 



B 



local ; 



B(r + 1)+B(r) 



(3) 



We calculate the structure functions of the local perpendicular mag- 
netic field component Bx, which corresponds to the Alfvenic fluc- 
tuations, at a variety of separations 1. 

The structure function values are binned according to scale 
parallel, h, and perpendicular, l±, to B] oca i. Nine linearly spaced 
bins are used in each direction covering the range 2,000 km to 
20,000 km, which is within, although towards the small scale end, 
of the inertial range. The result of this binning for one of the 65 
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Figure 2. second-order structure function of the perpendicular magnetic 
field component for one of the 65 solar wind intervals as a function of par- 
allel (iii) and perpendicular (/jj separation. 



intervals is shown in Fig. [2] It is representative of the average be- 
haviour, although in general each interval is more noisy and has 
less coverage than this. Most of the bin values in this figure are the 
average of a few thousand structure function values, although some 
(13%) are of a few hundred. It can be seen that the contours are 
elongated in the field parallel direction, indicating that the eddies 
are anisotropic with > feii . 

The data is also binned according to scale / and the angle 
Ob between 1 and Bi oca i. Nine linearly spaced scale bins are used 
over the range 2,000 km to 20,000 km and nine linearly spaced 
angular bins are used between 0° and 90°. Straight lines, in log- 
log space, are then fitted to the structure functions over the full 
scale range for each 6b bin and the power anisotropy is obtained 
by evaluating these fits at a scale of 10,000 km. The spectral in- 
dex in each 6b bin is found using the relation a = g+ 1, where 
— a is th e spectral index and g i s the structure function scaling 
exponent (Moni n & Yaglomll 19751) . This is similar to the work of 
lOsman & Horburvl (12009). except we bin the data with respect to 
the local field direction, since this appears to be the relevant mean 
field for the fluctuations, and we use many more intervals. 



2.3 Magnetic Field Anisotropy 

The results, averaged over all 65 slow wind intervals, are shown 
in Fig. [3] where the error bars are the standard error of the mean 
from averaging the intervals. They are similar to previous single 
spacecraft observations in the fast wind that show that power in- 
creases with 9b and that the spectral index varies fro m —2 at small 
angle s to between —5/3 and —3/2 at large angles l Horburv et al.l 
2008: iPodestH 12009: Luo & Wu1l2010l : Iwicks et al]|20ld 1201 ll) . 
This scale-dependent anisotropy, therefore, has now been seen in 
both fast and slow wind using two different measurement tech- 
niques. The power anisotropy is consistent with eddies elongated 
along the local magn etic field direction and wavevector anisotropy 



of the form k± > k\\ jChen et alfeOlObh . 

Th e —2 scaling at smal l Q B is cons istent with both the the- 
ories of lGoldreich&Sridhad i 19951) and iBoldvrevl §006), which 
describe cri tically balan ced Alfvenic turbulence. It has been sug- 
gested (e.g. Galtie jhoid) . however, that the parallel scaling of —2 
may be due to discontinuities in the data. Using the same technique 
at smaller scales in the dissipation range, the —2 scaling is not seen 
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Figure 3. Power anisotropy (upper) and spectral index anisotropy (lower) of 
the perpendicular magnetic field component in inertial range turbulence in 
the slow solar wind. The power anisotropy is calculated at / = 10,000 km. 
Spectral index values of —3/2, —5/3 and —2 are marked as dotted lines for 
reference. 



dChen et al hoiOal) because the physics of the turbulence is differ- 
ent at these scales. This suggests that the —2 scaling in the inertial 
range seen here is in fact due to the properties of the turbulence and 
not unrelated discontinuities. 

The perpendicular spectral index that we obtain here (for 
20° < 6 B < 90° ) is closer to -5/3 than - 3/2. This a grees with 
the pre diction of lGoldreich & Sridhaj (l995), rather than Boldvrevl 
(2006). Both of these theories, however, apply to balanced turbu- 
lence, i.e., cv = 0. As can be seen from Fig.Q] many of the inter- 
vals have large a c . This is common in the solar wind and various 
theories of imbalanced MHD turbulence have been proposed (e.g . 
i Lithwick et alj200llChandranll2008l:lBeresnvak & Lazariarj2o"o8l : 
IPerez & Boldvrevll2009l : |Podesta & BhattacharieekOld) , Differen- 
tiatin g between these, however, is beyond the scope of this paper 
Tsee lWicks et ai]j201 l|) for a recent observational test of these the- 
ories]. 

In the next section, we apply a similar analysis to RMHD sim- 
ulations. This enables a direct comparison to be made between tur- 
bulence in the solar wind and in numerical simulations. 



3 REDUCED MHD SIMULATIONS 
3.1 Simulation Description 

The RMHD equations, originally derived bv lStraussN 19761) . have 
been used previously to simulate various aspects of MHD turbu- 
lence (e.g. IPerez & Boldvrevl [20081 . 120091 : [Beresnvakil201lh . They 
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can be written in Elsasser potentials dSchekochihin et alj|2009h : 



Forced 



Decaying 



dt 



{c + ,vir}+{r,vic + }TV 2 ± {f+,r} 



(4) 



where {A, 6} = z ■ (Vj_A x V±B), i is the global mean field direc- 
tion, va is the Alfven speed and the Elsasser potentials are defined 
via 5z* = 8u± ± 8b± = z x Vj^*. 

Equations (O contain only the perpendicular fluctuations and 
are, therefore, suitable for simulating Alfvenic turbulence. They 
are also more efficient to simulate than MHD, since they involve 
only two scalar fields. Although originally derived from MHD, it 
has been shown that RMHD holds for a collisionless plasma such 
as the solar wind and may, the refore, be more generally appli- 
cable dSchekochihin et al] |2009). The RMHD derivation assumes 
anisotropy (k± S> hi ) and a strong mean field (So S> 5Bx), both of 
which are observed at the smallest scales of the solar wind inertial 
range. 

The simulation reported here solves the RMHD equations in 
a triply periodic cube of size (2k) 3 with a resolution of 512 3 . The 
Alfven speed is set to = 1 (making the Alfven crossing time 2n). 
It can be seen from equations l[4j that if the Alfven speed is scaled 
by a factor R and the z coordinate, which is the mean field direc- 
tion, is also scaled by R, the equations remain identical. This means 
that a given simulation corresponds to all values of R, and therefore 
all values of SB^/Bq if the box is also stretched in the z direction. 
The units of length in the perpendicular and parallel directions are 
independent of each other because the anisotropy is formally infi- 
nite and the fluctuation level is infinitely small under the RMHD 
asymptotic expansion. Different values of R can be chosen, setting 
the anisotropy and fluctuation level so that the same simulation can 
be compared to a variety of real world situations. 

The equations are solved pseudospectrally in x and y, and us- 
ing a centred finite difference scheme in z- The time step is chosen 
so that the Courant numbers based on both the fluctuation ampli- 
tude and the Alfven speed are much less than unity. With dissipa- 
tion and forcing terms the equations are 



1 



{r,vir}+{r,vi£ + }TVi{c + ,r} 



-vV5.(Vi£ ± ) + v l v2(viC ± )+.T 



(5) 



where v = 5 x 10~ 15 and v z = 1 x 10~ 4 are the viscosity coef- 
ficients and is the forcing term. In the x and y directions, a 
4th order hyperviscosity dissipation term is used, while in the z 
direction a very small Laplacian viscosity is added to prevent the 
high k z modes becoming unstable. Hyperviscosity is used so that 
the inertial range covers a wide enough range of scales to measure 
accurate scalings. The magnetic Prandtl number is Pr m = 1 and 
the initial conditions are a straight mean field with no fluctuations: 
b(r,r = 0) = x and u(r,f = 0) = 0. 

The simulation is initially forced on large scales (jfcj_ = 1,2 
and k z = 1) with Gaussian white noise forcing / , i.e. the random 
forcing amplitude is refreshed at each time step. This means that 
the input power can be controlled; it is set to unity in the code units 
to produce strong turbulence. We choose to force only the velocity 
to match possible sources of solar wind forcing, such as velocity 
shears or large scale Alfven waves, so / + = f~ at all times. We do 
not force the magnetic field since there is no known mechanism of 




Figure 4. Time series of simulation parameters: RMS variables, normalised 
cross helicity (J c , Alfven ratio and Fourier perpendicular spectral indices. 
For ( ^ 28 the simulation is forced and for f > 28 it is decaying. Spectral 
index values of —3/2 and —5/3 are marked as dotted lines for reference in 
the lower panel. 



breaking magnetic flux conservation at large scales. After a while, 
the forcing is removed and the simulation is left to freely decay. 

A time series of various simulation parameters is shown in 
Fig-El After the simulation begins, the values take a few time units 
to settle down, which is roughly the turnover time of the largest ed- 
dies. The transition between the forced and decaying periods of the 
simulation can be seen by the change in behaviour of all the quanti- 
ties at t = 28, marked by the dashed line. The top two panels show 
the root mean square (RMS) values of the Elsasser variables, veloc- 
ity and magnetic field. Their values up to t = 28 are determined by 
the forcing power and after / = 28 by the decay of the turbulence. 
One noticeable feature is the oscillation in the velocity and mag- 
netic field RMS values with a period ~ 2k. This is most lik ely due 
to large scale Alfven waves, also seen bv lBigotetaf] 12008), which 
should not significantly affect the average inertial range measure- 
ments. 

The normalised cross helicity <J C in the third panel is calcu- 
lated spectrally [equation (J)], as was done for the solar wind inter- 
vals in Section |2Tl and averaged over the range 7 ^ k^ 33. Dur- 
ing the forced period, o" c fluctuates above and below 0. When the 
forcing is re moved, \a c \ increases, as e xpected from dynamic align- 
ment theory dDobrowolnv e t al. 1980). The increase is fairly slow: 
\a c \ changes from 0.045 at t = 28 to 0. 13 at t = 78, which is consis- 
tent with previous decaying simulations (e .g. Grappi n et al. I I 19821 : 
iMatthaeus etail 19831 : iPouquet et all 19861 : lOughton et al.ll 19941) . 



The Alfven ratio r& is shown in the fourth panel, calculated 
over the same range as <J C . During the forced period 4 ^ t ^ 28, its 
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mean value is r& ~ 0.66, which is close to the solar wind observa- 
tions (TableQJ. As the turbulence decays, grows and approaches 
unity; this equipartitio n of energy is expected for MHD turbu- 
lence jKraichnanll 1 9651) . We note that the opp osite effect is seen in 
simulations without a str ong mean field (e.g. Ought on et al ] | 19941 ; 
iBiskamp & Mulled 19991) . in which the Alfven ratio decreases away 
from unity as the energy decays. The fact that the equipartition oc- 
curs only in the decaying period of our RMHD simulation, while 
solar wind observations show r& < 1 (e.g. [M atthaeus & Goldstein 
19821; iMarsch & Jul 1 199Cj ; IPodesta et al.ll2007l ; iBruno et alj|2007l ; 
Salem et al.l2009h . suggests that solar wind turbulence may be bet- 
ter described by a forced model. 

The perpendicular spectral indices for the velocity and mag- 
netic field are shown in the lower panel of Fig. [4] They are cal- 
culated from the gradients of the best fit lines to the perpen- 
dicular energy spectrum in log-log space over the range 7 ^ 
^ 33 every time unit. The perpendicular spectrum is calcu- 
lated as the sum of the energy in modes nearest to = ^Jk^ + ky 
for integer values of k±. During the forced period, the spec- 
tral indices are c loser to —3/2 than —5/3 , in agreement with 
previous results llMaron & Goldreich l fcoOll; iMtiller et al.l 120031; 



Miiller & Grappinl 120051; iMason et all 120081 ; IPerez & Boldvre 
2008; Grappin & Miiller 2010). When the forcing is removed, how- 
ever, they gradually steepen and appear to reach a steady value of 
—5/3 from t = 58 onwards. 

In the following analysis, we investigate the anisotropic scal- 
ing in the forced period 4 ^ t ^ 28, and the decaying period 
58 ^ t ^ 78. We assume that in each of these periods the turbu- 
lence is stationary and we can perform time averages over them. 
The averaged energy spectra are shown in Fig. [5] Before averaging, 
the decaying spectra are normalised so that the average energy over 
the range 7 ^ k^ ^ 33 for each is the same as that at r = 58. Gradi- 
ents of —5/3 and —3/2 are given for reference, although it is hard 
to tell the difference between these visually. It can be seen that for 
7 fex ^ 33 there are well defin ed power laws in a ll of the spectra. 
It has been suggested (e.g. lPerez & Bo ldvrev 2010) that the use of 
hyperviscosity may increase the bottleneck effect, altering the scal- 
ing. The spectra in Fig. [5] do not, however, display the increase of 
energy at small scales that is associated wit h the bottleneck effect 
and is seen in some MHD simulations (e.g. lCho & Vishniadl2000l ; 
lBeresnvakll201 lh . In the next section, we measure the anisotropic 
scaling using structure functions, which are expected t o be less sus- 
ceptib le to the bottleneck effect than Fourier spectra ( Dobl er et al.l 
120031) . 
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Figure 5. Perpendicular energy spectra of velocity (£""), magnetic field 
(E b ), and Elsasser variables (E ± ) in the forced and decaying periods of 
the simulation. Slopes of —5/3 and —3/2 are given for reference. 




Figure 6. second-order structure function of the perpendicular magnetic 
field component for one of the snapshots (t = 28) in the forced simulation 
as a function of parallel (iii) and perpendicular (/jj separation. 



3.2 Analysis Technique 

The technique we use to analyse the simulation data is similar to 
that used in Section |2~2l with modifications to account for the sim- 
ulation geometry. Firstly, the scaling factor R, which should be 
larger than unity for the RMHD equations to be valid, is chosen. 
Here, we set R = 4, which is a compromise between typical solar 
wind wavevector anisotropies k±/ku of between 2 and 3 and typi- 
cal 5B^/6o values of 0.1 (calculated from the data in Section l2~3l >. 
This means that the simulation, which was solved in a (2;r) 3 box, 
is now stretched to have a size (2jr) 2 x 8?r and the Alfven speed is 
set to 4. 

For a particular snapshot in time, many pairs of points in the 
simulation box are picked at random. The second-order structure 
function values of the local perpendicular velocity and magnetic 
field components are calculated and binned, as in Section l2~2l The 



structure function of the magnetic field binned with respect to In 
and Zj_ at t = 28 is shown in Fig. [6] There are on average 10 4 struc- 
ture function values in each bin. The structure function in Fig. [6] 
is representative of the general shape of the velocity and magnetic 
field structure functions in both the forced and decaying periods of 
the simulati on. Similarly to the so lar wind (Fig. |2j and previous 
simulations dCho & Vi shniac 200f3), the contours are elongated in 
the parallel direction. 

In the range 0.35 / 1.3, which corresponds approximately 
to 5 ^ k ^ 18, the structure functions are approximately power laws 
and we assume this to be the inertial range of the simulation. The 
spectral indices and the power anisotropy (calculated at / = 0.8) 
are found from the best fit lines to the data binned with respect to 
/ and 6b in this range. This is done for snapshots separated by 2 
time units, giving 13 snapshots for the forced period and 1 1 for the 
decaying period. 
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Figure 7. Power anisotropy (upper panels) and spectral index anisotropy (lower panels) of the velocity and magnetic field in the forced (left) and decaying 
(right) simulations. The power anisotropy is calculated at / = 0.8. Spectral index values of —3/2, —5/3 and —2 are marked as dotted lines for reference. 



3.3 Power and Spectral Index Anisotropy 

The power and spectral index anisotropy for the velocity and mag- 
netic field are shown in Fig. [7] The error bars are the standard error 
of the mean from averaging the results of the snapshots in each pe- 
riod. In both cases, the power increases with angle to the local mag- 
netic field 6b, as it does in the solar wind. For the forced case, the 
overall power in the magnetic field is larger than that in the veloc- 
ity, whereas in the decaying case they are similar. This is consistent 
with our previous discussion of the Alfven ratio being < 1 in the 
forced case and ~ 1 in the decaying case. It is also interesting 
to note that these curves are qualitatively similar in shape. No pre- 
diction for this shape has yet been made based on critical balance 
theory. 

In the forced case, there is a difference between the spectral in- 
dices of the velocity and magnetic field. The velocity spectral index 
varies from —2 at small B to —5/3 at 6b close to 90°. The mag- 
netic field spectral index is also —2 at small 6b but is less steep at 
larger 6b, having a value close to —3/2. The fact that both are steep 
at small angles shows that the turbulence is anisotropic and the 
—2 scaling is consistent wi th th e critical balance theories of both 
iGoldreich & Sridhaj J 19951) and iBoldvrevI d2006h . The difference 
at large angles, however, is unexpected, since theories of Alfvenic 
turbulence predict that both fields scale in the same way. 

For the decaying case, both fields show similar scaling. The 
spectral index is close to —5/3 for 6b close to 90° and much steeper 
at small 6 B : -2.33 ±0.03 for velocity and -2.30 ±0.03 for the 
magnetic field. Again, the steepening at low 6b shows the turbu- 
lence is anisotropic, although the spectra are steeper than the criti- 
cal balance prediction of —2. One possible explanation for the steep 
parallel scaling is that the turbulence may be transitioning to the 
weak regime, in which there is not thought to b e a parallel cascade 
dGoldreich & Sridharlll997l : lGaltier et~S1l2OO0l) .l Perez & Boldvrevl 



( 120081) observed the perpendicular spectral index steepening as the 
turbulence became weaker, and we may be observing a similar 
effect for the parallel index. In a different run (not shown here) 
that was forced less strongly, we observed overall steeper spectral 
indices at all angles. The perpendicular spectral index, however, 
seems to remain at —5/3 for many turnover times in the run here 
(Fig. rather than dropping to —2 as expected for weak turbu- 
lence. 

We now compare the spectral indices obtained through the 
structure function technique to the Fourier indices. The time series 
of the global perpendicular Fourier indices are shown in the lower 
panel of Fig. [4] It can be seen that both fields have spectral indices 
close to —3/2 during the forced period and then after a transition, 
reach a value of —5/3 in the decaying period. The mean values 
are -1.51 ±0.01 for the velocity and -1.47 ±0.01 for the mag- 
netic field in the forced period and —1.69 ±0.01 for the velocity 
and —1.653 ± 0.007 for the magnetic field in the decaying period. 
These are consistent with the perpendicular spectral indices mea- 
sured using structure functions, except for velocity in the forced 
period, which is close to —5/3. It is possible that this difference is 
caused by the forcing, which is localised at large scales in Fourier 
space, but may affect th e structure function, which m ixes small and 
large scale information (Davidso n & Pearsonll2005l) . 

The results we obtain here are broadly consistent with previ- 
ous simulations. Wavevector anisotropy of th e for m k \ > has 
been o bserved previously (e 
1 1994 iMatthaeus etd] Il99 



Shebalin et alJll983UOughton et al.1 



Milano et all 1200 if) . In particular. 



ICho & Vishniad (2000) observed a difference in anisotropic scal- 
ing between the velocity and magnetic field in their forced simu- 
lations. When the mean field was of a similar strength to the RMS 
fluctuations they obtained hi ~ k^ 1 for velocity but hi ~ k^ 5 for 
the magnetic field. 

We now compare the simulation and solar wind results. Firstly 
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we note that both sets of results are qualitatively similar. Power at a 
fixed scale is anisotropic and increases as 9b increases. All spectral 
index curves are anisotropic and steepen at small 9b as predicted 
by critical balance theories. The main difference between the solar 
wind and simulations is the value of the perpendicular spectral in- 
dex. For the magnetic field, we observe —5/3 in the solar wind and 
the decaying simulation but —3/2 in the forced simulation. Val- 
ues close to both —5/3 and —3/2 have been observed previously 
in the solar wind |Horburv et al1l2008l ; |Podestall2009l : iLuo & Wul 
l20ld : IWicks et alj 1201(1 1201 lh In both our forced and decay- 
ing simulations, the velocity has a perpendicular spectral index of 
—5/3. So lar wind measurement s , however, suggest that it is closer 
to -3/2 dMangenev et~ai]|200ll;|Podesta et all2007l ; iTessein et al.l 
120091 : ISalem et alj|2009l : IWicks et alj|201 lh . These differences in 
perpendicular spectral index remain an unsolved problem. 



4 LOCAL VS GLOBAL MEAN FIELD 

In this section, we investigate the difference between using the local 
and global mean magnetic field to define the parallel and perpen- 
dicular directions. Fig.[8]shows the spectral index anisotropy for the 
solar wind magnetic field in the upper panel, the forced simulation 
magnetic field in the centre panel and the forced simulation velocity 
in the lower panel. In each case, the results obtained using the local 
mean field are shown in green and those obtained using the global 
mean in orange. In the solar wind, the global mean field results are 
obtained by binning the structure function values according to their 
separation direction with respect to the average field of each inter- 
val. In the simulation, they are obtained by binning with respect to 
the average field over the whole simulation box (the z direction). 

In the solar wind, the results are similar for both local and 
global mean field methods, except for the error bar on the parallel 
scaling (0° < 9b < 10°). The error is larger in the global mean field 
case, which is partly due to the fact that there are fewer intervals 
where the global mean field is parallel to any separation vectors. 
This scaling, therefore, is less reliable and the error bar indicates 
that the data is marginally consistent with an isotropic spectral in- 
dex with respect to the global mean fie ld. To within errors, t hese re- 
sults are not inconsistent with those of lTessein et alj |2009j), which 
show that the spectral index is isotropic when measured with re- 
spect to the global mean field. 

In the forced simulation, it can be seen that when measured 
with respect to the global mean field, the spectral indices of both 
the magnetic field and velocity are much less anisotropic than when 
they are measured with respect to the local mean field, e.g. in 
the velocity at small 9b the spectral index is —1.76 ±0.02 using 
the global mean field compared to —1.97 ±0.02 using the local 
mean field. This is because the magnetic field fluctuations are large 
enough that the local mean field direction seen by an eddy is not 
the same as the global mean field direction. If the fluctuations are in 
critical balance, the angle between the local and global mean fields 
is / Bo ~ k\\ /kj_ . This suggests that when using the global mean 
field, the parallel scaling cannot be correctly distinguished from the 
perpendicular scaling, even for small 8B^/Bq, because the angle of 
measurement to the local mean field needs to be less than k«/kj_. 

This interpretation is in agreement with previous solar wind 
studies that have used local and global mean field methods. Those 
that use the global mean field method do not detect spectral index 
anisotropy dSari & Vallevlll976l : lTessein et alj|2009h and those that 
use a local mean field method do detect it ( Horburv et ai] |2008l : 
|Podestall2009l ; lLuo & Wu1l20ld:IWicks et al.ll2010L 1201 lb . A simi- 
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Figure 8. Spectral index anisotropy of magnetic field in the solar wind (up- 
per), magnetic field in the forced simulation (centre) and velocity in the 
forced simulation (lower) using the local and global mean field methods. 
Spectral index values of —3/2, —5/3 and —2 are marked as dotted lines for 
reference. 



lar situation is also seen in simulations, w here scaling anisotropy 
is detected when a local m ean field is used ( iCho & Vishniac 2000; 
iMaron & Goldreichl200ll) but not when a global mean field is used 
dGrappin & Mullej2010h . Here, we have shown that when keeping 
all other parameters constant, it is indeed the use of the global or 
local mean field that determines whether the anisotropic scaling is 
measured. It seems, therefore, that the Alfvenic fluctuations, both 
in solar wind turbulence and forced RMHD turbulence simulations, 
are more sensitive to the local mean field at the scale of the fluctu- 
ations than the global large scale field. 

In the decaying simulation (not shown in Fig.[8j», the local and 
global mean field methods are much more similar, with the paral- 
lel scaling being steeper than —2 in all cases. One possible reason 
for this is that the scale separation between the global mean field 
and the fluctuations is not large, meaning that the global and local 
mean fields are similar. This, combined with the smaller fluctua- 
tion amplitudes in the decaying simulation, could account for the 
observed behaviour. This could be tested by performing a decaying 
simulation with a larger inertial range. 



5 SUMMARY AND CONCLUSIONS 

In this paper, we measure the power and spectral index anisotropy 
of Alfvenic turbulence in the solar wind and RMHD simulations 
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using second-order structure functions. The analysis technique is 
essentially the same for both, allowing us to make a direct compar- 
ison. In the slow solar wind, we find that the magnetic field power 
and spectral index are anisotropic with respect to the local mag- 
netic field direction. This anisotropy has now been seen by several 
different methods in both fast and slow wind. In both forced and de- 
caying simulations we also find that the power and spectral index 
are anisotropic in both the velocity and magnetic field. 

In the solar wind, the perpendicular spectral index of the 
magnetic field is close to — 5/3, in agreement with the theory of 
Goldr eich & SridhaJ Jl995h . In the forced simulation, the perpen- 
dicular spectral indices are close to —5/3 for velocity and —3/2 
for the magnetic field. We are not aware of any theory that can ac- 
count for this difference, although it may be caused by the velocity 
forcing. In the decaying simulation, the perpendicular spectral in- 
dex is close to —5/3 for both the velocity and magnetic field. In all 
cases, the spectral index steepens at small angles to the magnetic 
field. The parallel scaling obtained in the solar wind and forced 
simulations is close to —2, which agrees with the the orie s based on 
critica l balance of both iGoldreich & Sridhar] d 19951) and lBoldvrevI 
(2006). The parallel spectral indices in the decaying simulation are 
-2.33 ±0.03 for the velocity and -2.30 ±0.03 for the magnetic 
field, which are steeper than the critical balance predictions. 

We also find that when measuring the anisotropy of the fluctu- 
ations in the forced simulation with respect to the global magnetic 
field, rather than the local mean field, the spectral indices are much 
less anisotropic. This is expected for critically balanced turbulence 
and is also consistent with previous solar wind and simulation re- 
sults: those that used the local mean field saw anisotropic scaling 
and those that used the global mean field did not. 
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